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Abstract. It is well established now that the solar atmosphere, from photosphere to the corona and the 
solar wind is a highly structured medium. Satellite observations have confirmed the presence of steady 
flows. Here, we investigate the parallel propagation of magnetohydrodynamic (MHD) surface waves trav- 
elling along an ideal incompressible flowing plasma slab surrounded by flowing plasma environment in the 
framework of the Hall magnetohydrodynamics. The propagation properties of the waves are studied in a 
reference frame moving with the mass flow outside the slab. In general, flows change the waves' phase 
velocities compared to their magnitudes in a static MHD plasma slab and the Hall effect limits the range 
of waves' propagation. On the other hand, when the relative Alfvenic Mach number is negative, the flow 
extends the waves propagation range beyond that limit (owing to the Hall effect) and can cause the trig- 
gering of the Kelvin-Helmholtz instability whose onset begins at specific critical wave numbers. It turns 
out that the interval of Alfvenic Mach numbers for which the surface modes are unstable critically depends 
on the ratio between mass densities outside and inside the flux tube. 

PACS. 96.50.Tf MHD waves; plasma waves, turbulence - 96.50.Ci Solar wind plasma; sources of solar 
wind 
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1 Introduction 

Various waves and oscillations which occur in structured solar atmosphere were intensively studied over the past 
three decades and an exclusive review of their theory the reader can find in [l] and references therein. Next step in 
studying the wave phenomena in solar and stellar atmospheres was the consideration of steady flows there. Satellite 
measurements of plasma characteristics of, for instance, the solar wind and coronal plumes flows, such as the magnetic 
field, the thermal and fiow velocity and density of plasma or plasma compositions, arc important to understand the 
various plasma wave modes which may arise. However, wave analysis requires further information and special tools to 
identify which set of modes is contributing to observed wave features. In practice, one may use filters to perform the 
so-called pattern recognition to detect the various kind of modes that may propagate in plasma and to determine their 
contribution to the wave energy [2]. Another important issue is the question for the waves' stability. The magnetosonic 
waves in structured atmospheres with steady flows have been examined by Nakariakov and Roberts [5], Nakariakov 
et al. [1], Andrics and Goossens [S], and Terra-Homem et al. [B]. Andrics and Goosscns studied also the conditions 
at which resonant flow and Kelvin-Hclmholtz instabilities take place, and Terra-Homem et al. have investigated the 
effect of a steady flow on the linear and nonlinear wave propagation in various solar structures. 

It is worth mentioning that all the aforementioned studies were performed in the framework of standard magne- 
tohydrodynamics. It was Lighthill |7| who pointed out almost 50 years ago that for an adequate description of wave 
phenomena in fusion and astrophysical plasmas one has to include the Hall term, TOi(j x B)/e/9, in the generalized 
Ohm's law. That approach is termed Hall magnctohydrodynamics (Hall MHD). In this way, it is possible to describe 
waves with frequencies up to w ujci- Because the model still neglects the electron mass, it is limited to frequencies well 
below the lower hybrid frequency: uj <C wm- Generally speaking, the theory of Hall MHD is relevant to plasma dynam- 
ics occurring on length scales shorter than an ion inertial length, L < Zuaii ~ c/wpi (where c is the speed of light and 
ojpi is the ion plasma frequency), and time scales of the order or shorter than the ion cyclotron period (t < utd^^) [8]. 
Thus the Hall MHD should affect the dispersion properties of MHD waves in spatially bounded magnetized plasmas. 

The first paper devoted to the propagation of fast MHD and ion-cyclotron surface waves at a plasma-vacuum 
interface in the limit of cold plasma was that by Cramer and Donnelly |^. Later on, Cramer |10| generalizes that model 
obtaining the dispersion of nonlinear surface Alfvcn waves. It is worth noticing that even in unbounded magnctoplasma 
at purely parallel propagation the Alfven waves become dispersive when the Hall current is taken into account in the 
basic MHD equations [TI]. In Ref. [T^] the Hall MHD has been applied in studying the parallel propagation of fast 
wave in an ideal static plasma slab. There were derived four boundary conditions (see Sec. HI in |12j ) necessary for 
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obtaining the dispersion relations of sausage and kink modes in spatially bounded magnetized plasmas. That approach 
has been applied in studying the oblique propagation in the same geometry both for low-/? and finite-/? plasmas |13|14j . 

The first study on surface- wave parallel propagation in a flowing ideal MHD flux tube surrounded by a static plasma 
environment (both embedded in a constant magnetic field Bg) in the framework of the Hall MHD was performed by 
Miteva et al. |15| . It has been shown that while in a static plasma slab the hydromagnetic surface waves (sausage and 
kink modes) are Alfvcn ones (their phase velocities are close to the Alfvcn speed in the layer), in slabs with steady 
flows they become super-Alfvenic waves. Moreover, as it is logical to expect, there exist two type of waves: forward 
and backward propagating ones, bearing in mind that the flow velocity deflnes the positive (forward) direction. An 
in-deep examination of wave modes in flowing solar-wind-flux-tube magnetized plasmas for finite-/? and zero-/? ionized 
media has been performed in |16ll7j (see also the references in those papers). 

Hall MHD is relevant not only to linear MHD waves, but also to nonlinear ones |18I19I20I21| . Dispersive effects 
caused by Hall currents perpendicular to the ambient magnetic field can influence the generation and propagation of 
shock waves Recently Clack and Ballai [53] studied the nonlinear wave propagation of resonant slow magnetoa- 
coustic waves in plasmas with strongly anisotropic viscosity and thermal conductivity alongside of dispersive effects 
due to Hall currents. They show that the nonlinear governing equation describing the dynamics of nonlinear resonant 
slow waves is supplemented by a term which describes nonlinear dispersion and is of the same order of magnitude as 
nonlinearity and dissipation. In the case of stationary nonlinear waves the Hall-MHD equations can be rewritten in 
the so-called Sakai-Sonnerup set of equations that describe the plasma state and provide oscillatory and solitary types 
of solutions ^24] . The overall parameter study on the polarization characteristics, together with the magnetic field 
components and density variations in different ranges of solutions performed in that paper might be further on applied 
to the theoretical treatment of particle interaction with such waves, e.g., at shocks in space plasmas, possibly leading 
to particle acceleration. MHD parametric instabilities of parallel propagating incoherent Alfven waves are influenced 
by the Hall effect and that is especially important for the left-hand polarized Alfven waves [25] . Ruderman and Caillol 
pS] claim that the left-hand polarized Alfvcn waves in a Hall plasma are actually subject to three different types of 
instabilities, namely modulational, decay, and beat instability, while the right-hand polarized waves are subject only to 
decay instability. A new trend in the Hall MHD is its application to partially ionized plasmas j^T] - the Hall effect may 
play an important role in the dynamics of weakly ionized systems such as the Earth's ionosphere and protoplanetary 
discs. The exact nonlinear cylindrical solution for incompressible Hall-MHD waves, including dissipation, essentially 
from electron-neutral collisions, obtained in a uniformly rotating, weakly ionized plasma, derived by Krishan and 
Varghese |28) . demonstrates the dispersive nature of the waves, introduced by the Hall effect, at large axial and radial 
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wave numbers. Such short-wavelength and arbitrarily large-amplitude waves could contribute toward the heating of 
the solar atmosphere. One of the most important mechanisms for solar atmosphere's heating and especially the solar 
wind is the wave turbulence, which is also affected by the Hall currents [29130131132133] . The strongest competitor of 
the wave turbulence heating mechanism, notably the magnetic rcconncction, is similarly heavily influenced by the Hall 
effect |34l35l36l37l38j . One even speaks for Hall reconnection being akin to Petschck rcconncction model. It is worth 
mentioning that the Hall effect can influence the magnetic field dynamics in dense molecular clouds [39] , magnetoro- 
tational instability [40], star formation [41], as well as compact objects [42]. As seen, the Hall MHD has impact on 
many important astrophysical phenomena and objects. 

Here, we investigate the influence of flow velocities on the dispersion characteristics and stability of hydromagnetic 
surface waves (sausage and kink modes) travelling along an infinitely conducting, magnetized jet moving past also 
(with a different speed) infinitely conducting, magnetized plasmas. If in the solar corona plasma /3 (the ratio of gas to 
magnetic pressure) is much less than unity, in the solar wind flux tubes it is /3 ~ 1. Since we are going to study the 
wave propagation in flowing solar wind plasma, we can assume that we have a 'high-/3' magnetized plasma and treat it 
as an incompressible fluid. For simplicity, we consider a planar jet of width 2x0 (embedded together with environments 
in a constant magnetic field Bq directed along the z axis), allowing for different plasma densities within and outside 
the jet. Pa and pc, respectively. The most natural discontinuity which occurs at the surfaces bounding the layer is 
the tangential one because it is the discontinuity that ensures a static pressure balance |43| . For typical values of the 
ambient constant magnetic field _Bo = 5 x 10^^ T and the electron number density inside the jet rio = 2.43 x 10^ m""^ 
at 1 AU |43| . the ion cyclotron frequency tJci/STr = 76 mHz, and the Hall scale length (= WAo/^ci, which is equivalent 
to c/wpi) is ^Haii ~ 150 km. This scale length is small, but not negligible compared to layer's width of a few hundred 
kilometers. Here, we introduce a scale parameter e — ^Haii/a^o called the Roll parameter. In the limit of £ 0, the Hall 
MHD system reduces to the conventional MHD system. Our choice for that parameter is £ = 0.4. The flow speeds of 
the jet and its environment are generally rather irregular. For investigating the stability of the travelling MHD waves 
it is convenient to consider the wave propagation in a frame of reference attached to the flowing environment. Thus 
we can define the relative flow velocity U''''' = Uq — Uc (C/q and [4 being the steady flow speeds inside and outside 
the flux tube, respectively) as an entry parameter whose value determines the stability/instability status of the jet. 
As usual, we normalize that relative flow velocity with respect to the Alfvcn speed in the jet and call it Alfvcnic 
Mach number Ma, omitting for simplicity the superscript "rel". Another important entry parameter of the problem is 
r] ~ Pc/ Po- It turns out that the waves' dispersion characteristics and their stability critically depend on the magnitude 
of rj. Our preliminary numerical studies, based on a criterion for rising the Kelvin-Helmholtz instability (namely that 
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the modulus of the steady flow velocity must be larger than the sum of Alfvcn speeds inside and outside the flux tube) 
proposed by Andries and Goossens [5], showed that in order to expect instability onset at some reasonable values of 
the relative Alfvcnic Mach number Ma, one should assume fairly large values of the parameter 77. In this study we take 
77 = 10 which means that the Alfvcn speed inside the jet is roughly three times larger than that in the environment - 
actually an "edge" choice, but still acceptable. Thus the waves' dispersion characteristics (the dependence of the wave 
phase velocity Wph = tj/fc on the wave number k) and their stability states are determined by the three parameters rj, 
e, and A/a, two of which are flxed (77 and s) and the third one. Ma, is running. 

The basic equations and wave dispersion relations will be exposed in Sec. 2 of the paper, and the numerical results 
and discussions - in Sec. 3. Section 4 summarizes the new flndings and comments on future improvements of this 
study. 



2 Basic equations and dispersion relations 

The jet's interfaces are the surfaces x — ±xq, the uniform magnetic field Bg and the steady flow velocities Vq.c point 
in the z direction. The wave vector k lies also along the z axis and its direction is the same as that of Bq and Uo,e- 
As we have already mentioned we will study the waves propagation in a reference frame affixed to the environment. 
Thus the steady fiow velocity in the slab is U = U''''' = Uq — Uc and zero outside. The basic equations for the 
incompressible Hall-MHD waves are the linearized equations governing the evolution of the perturbed fluid velocity 
Sv and wave magnetic field (5B: 

p^5v + p(U ■ V)(5v + V(— Bo • SB) - — (Bq • V)SB = 0, (1) 
at fiQ fio 

a 2 

— 6B + (U • V)(5B - (Bo • V)<5v + BoV • <5v + ^z • VV x <5B = 0, (2) 

Ot UJci 

with the constrains 

V • = 0, (3) 

V • (5B = 0, (4) 

where va = So/(mo/o)^^^ is the Alfvcn speed and ^0 is the permeability of free space. After Fourier transforming all 
perturbed quantities ozg{x) exp(— icjt + ifcz), we derive two coupled second order differential equations for 5vx and 5vy, 
notably 



d^ , .a -I 



fc2 _ i k^^^^ ^ 
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and 



where 



e = , and a ~ 



Wei V ^^A 

with Wei = P-Bo/m\ being the ion angular cyclotron frequency (rrii is ion/proton mass and e is the elementary electric 
charge). Note that we have different es and as inside and outside the jet. We seek the solutions to coupled equations 
([5]) and ([5]) in the form 

5v.j,[x) = f [exp(-Ka;) ^ exp(K.T)] , 
Svy{x) = ih [exp(— Kx) =F exp(Ka;)] , 
anticipating surface waves with attenuation coefficient k, and obtain the following set of equations 

(n^ - f - - k^) h ^ 0. 

This set of equations yields the following expressions for k: 

Ki = k, 

2 / 21 1/2 



K2 



fc [1 - (a - 1) = 



That means that there are in fact two pairs of attenuation coefficients: (fc, TOq) mside the flux tube and (fc, TOc) outside 
the jet, respectively. 

As is known, on a bounded MHD plasma wave guide (cylinder or layer) two types of waves may exist. Recall that for 
a slab geometry the general solutions to the equations governing 6vx and Svy are sought in the form of superpositions 
of cosh KqX and sinh KqX. Those solutions contain waves whose shape is defined by the cosh function (they are called 
kink waves) and another type of waves associated with the sinh function (which are termed sausage waves). The 
transverse structure of both waves inside the slab is determined by the two attenuation coefficients Koi,2 (i.e., k and 
mo). Thus the solutions for Sv^ and Svy inside the slab (\x\ < xq) assuming a sausage wave form, accordingly, are 

sinh kx sinh nioX 

ovx[x) = /i . , , h h- 



sinh kxo sinh moXo 
and 

„ . . . J, „ sinh kx sinh nioX 

OVy{x) ^ i/iGoi . , , 1-1/2^02^-; , 

smh kxQ smh rrioXQ 
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where 

Q — ° ol,2 

These solutions have been obtained by imposing the compatibiUty condition for the set of equations for / and h as 
unknown quantities. We also notice that Goi = and Go2 = (lo — l)/co- For a kink surface- wave form the expressions 
for perturbed fluid velocity components have the same description - it is only necessary to replace sinh with cosh. The 
solutions outside the layer (identical for both modes) are 



and 



SVy{x) 



Here, as above, 



ai exp [—k{x — xq)] + a2 exp [—mc{x — xq)] for x > xq, 
Pi exp [k{x + xq)] + [i2 exp [inc{x + xq)] for x < — .tq, 

iaiGoi exp [—k{x — xq)] + ia2Gc2 exp [— 77io(x — xq)] for x > xq, 
i/?iGoi exp [k{x + xq)] + i/32Go2 exp [toc(x + xq)] for x < —xq. 



^ _ Co ^cl.2 



flc — 1 k"^ ' 

and, similarly, Goi = and Gc2 — [o-c — l)/co- 

Having derived the expressions for the perturbed fluid velocity components Svx and Svy, one can calculate the 
perturbed total pressure, which in our case is the perturbed magnetic pressure only, and arrive at 

Sptotai{x) = —Bo • SB = i ^^——v\ 



liQ w — k • U 



(u; - k ■ U) 

k'^v\ 



2 



d , , , . w-k-U d 

—dvx{x) + 1 — 

da; (jJci ax 



The perturbed wave electric field can be obtained from the generalized Ohm's law 



mi 



E = -V X B H X B, 

ep 

which, by means of Ampere's law (multiplied vectorially by B) yields 



(5E = Bo X ( 5v - /Haii^V X (5B ) - U X (5B. 



Its components are 

^F,^(x] = - 

a; — k • U V kvA 



SE.ix) = —^(^^-^yBoSu,ix), 
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and 



5Ey{x) 



w-k-U 

The above expressions for perturbed quantities are used when implementing the boundary conditions. 

It follows from the solutions to the basic equations ^ and for the perturbed fluid velocity components Svx and 
Svy that the number of integration constants is six. However, because of the symmetry (or antisymmetry), the two 
/?i,2 amplitudes are in fact directly obtainable from the x > xo solutions ~ indexed as and f3s are not independent. 
Thus we can derive the dispersion relations by applying only four boundary conditions at one interface, for example, 
at X = Xo- These boundary conditions, as we already mentioned in Sec. (TJ are derived and discussed in \12\ . We 
can borrow them except the first one, the continuity of Svx across the interface, as in the present case of a flowing 
slab this condition must be replaced by the continuity of Svx/ {uj — \i ■ V) [H]. The rest of the boundary conditions 
are the continuity of the perturbed pressure (5ptotai, the j/-component of perturbed wave electric field 6Ey, and the 
x-component of perturbed electric displacement 6Dx = Sq {K^x^Ex + KxySEy) (where Sq is the permittivity of the 
free space) at the interface. In the last boundary condition, Kxx and Kxy are the low-frequency components of the 
plasma dielectric tensor [45] 



K. 



vi 



and Kx 



.W-k-U c2 



By imposing the boundary conditions at the interface x = Xq, and after some straightforward algebra, we finally 
arrive at the dispersion relations for parallel propagation of sausage and kink waves in a planar jet, surrounded by 
steady plasma media 



k-U 



kVAo 



- 1 



( ^ 



Pa \kvAo 



- 1 



( tanh\ 
icoth)^"° 



(7) 



1 + 0) 



~ 2 Pe / tanh\ 



Oo vcoth j 



kXr 



1 - UJPc/Po 
1-Q {Pc/Pof 



where 



and 



k U 



CJ — k • U LUci 

As can be seen, the wave frequency lu is Doppler-shifted inside the jet. These dispersion relations can be extracted 
from Eq. (12) in [16] in the limit Cg — ^ oo (sound speed much larger than the Alfvcn one) by letting Uc = and 
Uo = U. When the flux tube is static (U = 0), Eq. ([7]) coincides with Eq. (18) in Ref. [12]. Finally, when one neglects 
the Hall effect (eo = 0) one gets the well-known, derived by Edwin and Roberts [17], dispersion relations in the limit 
Cs — > OO. We should also emphasize that in an incompressible jet both modes are pure surface waves. A question which 
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immediately arises is what type of MHD waves arc those described by Eq. ([7|)? From the three well-known MHD 
linear modes propagating in infinite compressible magnetized plasmas (Alfven, fast and slow magnetosonic waves) in 
incompressible limit survive only two, notably the shear and pscudo Alfven waves. The latter is the incompressible 
vestige of the slow mode of compressible MHD. The displacement vector of a shear Alfven wave is perpendicular to 
the plane defined by its wave vector, k, and a uniform background magnetic field, Bq, whereas that of a pseudo Alfven 
wave lies in this plane. The two wave modes share the dispersion relation 



LiJ 



2 vk-Bo)2 _ , s2 

- ■ (K||WAj 



and propagate with group velocity, va, either parallel or antiparallel to Bq depending upon the sign of It seems 
that our Hall-MHD modes travelling along the jet are akin to the pseudo Alfven waves. 
The dispersion relation of each mode can be symbolically written down in the form: 

fc, parameters) = 0, (8) 

where the function argument 'parameters' includes data specific for the jet - they will be listed shortly. As we are 
interested in the stability of the surface waves running at the jet interfaces, we have to assume that the wave frequency 
is complex, i.e., w ^ cj + 17, where 7 is the expected instability growth rate. Thus dispersion equations ([7]) become 
complex and their solving is not a trivial problem [3^. When studying dispersion characteristics of MHD waves, one 
usually plots the dependence of the wave phase velocity Wph as function of the wave number fc. For numerical solving 
of equations ([7]) we normalize all quantities by defining the dimcnsionless wave phase velocity = w/fcuAo, wave 
number K ~ kxQ, and the (relative) Alfvcnic Mach number A/a ~ U/vao, respectively. Accordingly 

[{u-k- U)/(fci;Ao)]' = (l^ph - Ma)' , c^V(fc«Ao)' - V^L 

oj = VpYi/ {Vph- Ma) and ^ K {Vph - Ma) l-Raii/xo- 

Note, that /naii/a^'o = £ (alongside with rj and A/a) is an entry parameter which has to be specified at the start of the 
numerical procedure. Thus we have to solve normalized dispersion relations ([7]), having now the form 



(F,H-M.f-l + („K^.-l)(^-M 



K 



(9) 



2^2 



K 



l/ph(l - r;) - A/a _p 



yph(l - r?2) - A/a 

Recall that we consider the normalized wave phase velocity Vph as a complex number and wc shall look for the 
dependencies of the real and imaginary parts of V^h as functions of the real dimcnsionless wave number K at given 
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values of the three entry parameters. It can be easily seen from above equations that they are cubic ones with respect 
to V^h- Hence, the dimensionless dispersion relation of, for example, the kink mode can be displayed in the form: 



AVX + SVpi + CFph + Z? = 0, 



(10) 



where 



A = {l + T]cothK){l - T]){l + T]^ e^K^), 

B = -Ma [2(1 - T]){1 + T]- e^K^) + (1 + ?7Cothii:)(l - e^K^)] , 
C = M'l [(1 - r/)(l + ?7 - e^K^) + 1 - e^K^] - (1 + coth A')(l - r]^), 
D = Ma [~MI{1 - e^K^) + 1 + coth A'] . 

The dimensionless dispersion relation for sausage mode is similar - one has simply to replace coth by tanh. 
Cubic equations will be solved on using Cardano's formulas [48]. First let us define a variable /: 



Next we define g: 



Finally we define h: 



_ 2 BC D 

^ ^ 27^ ^ A' 



If /i > 0, there is only one real root and two complex conjugate ones. When /i ^ 0, all three roots are real. 
When h < 0, the real solutions to the cubic Eq. PH)) are: 

Vi =2\ cos -, 

V 3 3 3A' 



(11) 



1^2,3 = -2a/- 



— cos — ± — , 

3 V3 3/ 3A' 



(12) 
(13) 



where 



2V^(773F. 

For simplicity, we have dropped the subscript 'ph' to the normalized phase velocity V . 
When /i > 0, the real root of Eq. pO]) is given by 



Vq = M + N 



B 

3l' 



(14) 
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where 

+ and N ^ i^-g /2 - -/h^'^ . 



M 

The two complex roots are given by ± where 



V, = -\{M + N)~^, (15) 



and 



V,^^{M-N). (16) 

During the numerical solving of dispersion equations Q we can look what sign possesses h, defined by Eq. (jlip . 
and in the ranges of the normalized wave numbers K, where h is positive, we can expect complex solutions, i.e., 
amplification or damping of the waves due to their interaction with the flow. 



3 Numerical results and discussion 

Before starting the numerical solving of dispersion equations © we have to specify the entry jet's parameters [c.f., 
Eq. ([5])]. As was told in the Introduction section, we take the Hall parameter e ~ 0.4. The relative Alfvcnic Mach 
number, Ma, will be a running parameter, from (for a static flux tube) to some reasonable values. These values 
to some extent depend upon the choice of the third parameter, 77, equal to the ratio of plasma densities outside and 
inside the jet. When studying Kelvin -Hclmholtz instabilities of MHD waves in the coronal plume-interplume region 
in the framework of standard magnetohydrodynamics, Andries and Goossens [5] find that in the /3 = case one can 
expect that the instability will occur approximately for \U\ > uao + VAo- As we have already mentioned in Sec.[Tl we 
assume that this estimation is valid for our case of an incompressible plasma jet. After normalizing all velocities with 
respect to the Alfven speed inside the slab, uao, we get 

VAo VAo 

or 

|Ma|>1 + ^. (17) 

As seen from above inequality, for some small rjs the relative Alfvcnic Mach number, AIa, might become rather large 
in order to register the onset of Kelvin-Hclmholtz instability. Such large Mach numbers imply that for magnitudes of 
the Alfven speed inside a solar wind jet in the range, say, of 60-100 kms~^ the difference between Alfven speeds in 
the jet and its environment should be of the order of a few hundred kilometers per second which is unlikely to occur. 
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Moreover, for relatively large Alfvcnic Mach numbers dispersion curves of Hall-MHD surface modes become rather 
complicated. That is why we choose 77 = 10, which means that one can expect a onset of instability at \Ma\ > 1.316 
that looks as a reasonable value. 

A specific feature of the surface Hall-MHD waves travelling along an incompressible static plasma layer is that 
there exists a limiting dimensionless wave number Kn^n beyond which the wave propagation is no longer possible. 
That limiting wave number is given by [12j 

Knn.it = il + v)'^V£- (18) 

For our choice of rj and e, -R'limit ~ 8.292. With approaching that wave number the wave phase velocity becomes very 
large. It is interesting to sec whether the steady flow will change that limiting value. 

Let us first start with the kink mode. For solving the corresponding dispersion relation ^ we generally use the roots 
of the cubic equation given by Eqs. (fT^ - P^ . We begin the numerical solving with the relative Ma = (corresponding 
to a static fiux tube) running the dimensionless wave number K from 0.005 to 10. As naturally to expect, in that wave 
number region h is negative and we get real values for the normalized phase velocity - the corresponding dispersion 
curve is labeled by '0' in Fig. [TJ As seen from that figure, the kink wave is generally a sub- Alfvcnic one; at K w 7.2 
its phase velocity becomes equal to the Alfven speed and starts quickly to grow up reaching rather large values (up to 
800 even more) whence K Kn^in- What is going on as the layer possesses any flow velocity? If we take, for instance, 
Ma = 0.5, one can construct from the roots of the cubic equation two dispersion curves, both being real solutions, 
however, one curve with a positive phase velocity, and another curve with a negative one. The first dispersion curve 
as seen from Fig. [5] is slightly above the 'O'th curve, while the second dispersion curve starts with the negative value 
of —0.32 which becomes large in magnitude at K approaching its critical value, and what is more interesting, it goes 
beyond the i^nmit, now decreasing in magnitude. In this extended propagation range there exists a complex solution 
to the cubic equation with positive imaginary part, i.e., the wave becomes unstable. However, bearing in mind that it 
is unlikely to observe/detect such a backward wave in a solar wind tube (supposing that the wave phase and group 
velocities have the same direction), we have to drop all the solutions with negative phase velocities. The dispersion 
curves of the second type will be discussed in another, physically acceptable situation, soon. 

With increasing the relative Alfvcnic Mach number the dispersion curves with positive velocities initially lie below 
the 'O'th curve, but afterwards, for values of the normalized wave number between 3 and 5.5, they cross the curve 
with label '0' staying on the left side of that curve. In other words, the flow velocity slightly diminishes the value of 
the -fCiimit without allowing them to propagate beyond those limiting A's. Note also that the course of the dispersion 
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curves is not monotonous - initially, for small A/as, they lie above the neutral ('O'th labeled) curve in the range of 
small and average dimensionless wave numbers, while for Ma ^ 1 those curves set up below neutral dispersion curve. 

For any negative value of the relative Alfvenic Mach number Ma we have as before two set of solutions. The 
real ones are negative and represent a mirror image of dispersion curves shown in Fig. [1] We do not plot them for 
the same reason; although mathematically correct they are not acceptable from a physical point of view. The most 
interesting case are the dispersion curves shown in Figs. [3] and H) It is clearly seen that for each Ma, in fact, two 
distinctive dispersion curves merge at K « A'umit- The solutions to the dispersion relation for the curves lying on the 
right side of the merging vertical line at some A"s become complex with positive imaginary part, i.e., there the waves 
are unstable. The other important observation is that now the waves' phase velocities do not grow too much when 
K A'liiiiit. That is especially true for the dispersion curves associated with relatively large in modulus A/as - see, for 
example, the dispersion curve labeled by '—1.45.' Another important feature is the circumstance that all dispersion 
curves lying on the left side of the merging vertical line correspond to a stable (generally with complicated shapes of 
the dispersion curves) waves' propagation! We also note that parts of the dispersion curves, corresponding to a stable 
wave propagation, continue smoothly beyond the /•sTiimit ~ see, for example, in Fig. |4] the dispersion curve labeled by 
' — 1.25' that ends at kx^ = 10 with Wph/^Ao = 0.13. It turns out that the Hall current makes the waves stable in 
the wave number propagation range between and ifumit an instability onset is only possible at some critical wave 
numbers larger than A'limit- The growth rates of kink waves in the instability region are plotted in Fig. [Sj We would 
like to emphasize that all complex solutions were checked for some selected dimensionless wave numbers on using the 
complex versions of Newton-Raphson and Miiller [IS] methods. 

The case with sausage Hall-MHD waves is similar in many ways, although there are some specific features. First, 
for positive relative Alfvenic Mach numbers the real solutions in the long-wavelength limit (small Ks) are much more 
complicated. That can be seen in Figs. [B] and [T] The most striking issue is the existence of a bordering dispersion curve 
(with label '1.315') which divides the rest dispersion curves into two types. The dispersion curves with A'/a < 1.315 
begin as super- Alfvenic waves with decreasing phase velocities which passing through a minimum start to grow and at 
around K « 5 cross the neutral dispersion curve (that corresponding to a static slab). After that they quickly increase 
their magnitudes reaching very large values at K —f ifiimit • The second type of dispersion curves consist of two families 
of curves: ones at very narrow regions of small A's, and others starting with negligibly small negative phase velocities 
which later on passing through inverted "s-shaped" parts continuously increase their speeds reaching large values at 
K — !■ ATiiiiiit. It is worth noticing that the second type's dispersion curves possess multiple values at a fixed K in the 
range between and 1.22. 
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The most intriguing question is how the dispersion curves will behave as A/a is negative? The answer is illustrated 
in Figs. [8] and m ActuaUy there is no a big surprise - more or less the dispersion curves are similar to those of the kink 
mode. However, we should immediately notice two differences: (i) the dramatic changes in the shapes of the dispersion 
curves start at Mp^ = —1.3 (vs. —1.25 for the kink mode), and (ii) parts of dispersion curves for A/a ^ —1.3 with 
dimensionless wave number between and 8.2 are negative. All these forward and backward sausage waves are stable. 
Unstable are only those waves (like for the case of kink waves) whose dispersion curves are on the right of the merging 
vertical line (look at Fig. [8|) . The growth rates of such unstable sausage Hall-MHD modes are shown in Fig. [10] and 
each of them starts at some critical normalized wave number. 

We should recall that our theory is a linear one and since at the instability onset the waves amplitudes begin rising, 
for a further waves' evolution one must employ a nonlinear approach. Nevertheless, results, obtained here, can be used 
as a start point for a deeper investigation of the wave propagation in flowing solar flux-tube plasmas in the framework 
of the Hall magnetohydrodynamic. 

An interesting issue which springs to mind is how the instability growth rate depends on the value of the entry 
parameter 77, say, at a fixed dimensionless wave number. Let us do such an examination for the kink mode and let 
our choice for the fixed wave number be (see Fig. [S]) kxo = 8.5. The wave growth rate depends not only on K but 
also on the relative Alfvenic Mach number Ma. In Fig. [11] we show a family of curves depicting the dependence of the 
normalized wave growth rate as a function of the plasma densities ratio rj ~ Pel Po for various Alfvenic Mach numbers 
between —1.45 and —0.5. As seen, one observes two local maxima: one aX rj = \ and another around rj ~ 10. The first 
maximum is not surprising - it is well known that the Kelvin-Helmlioltz instability is easily excited when the densities 
of the two adjacent flowing media are approximately the same. The only exception here is the curve corresponding 
to A/a = —0.5 - a maximum of the wave growth rate for that value of the Alfvenic Mach number one can expect 
for values of bigger than 10. The second local maxima are obviously depending on the magnitude of the relative 
Alfvenic Mach number A/a. 

Another curious question is how the propagation and stability properties of the Hall-MHD surface waves change 
with the value of 77. If we take 77 = 4, the dispersion curves of kink waves for negative values of the relative Alfvenic 
Mach number are plotted in Fig. 1121 As seen, the picture is similar to that shown in Fig.]!] however, with an distinctive 
feature, notably the curves corresponding to A/a < — 1 represent stable wave propagation even on the right side of the 
vertical merging line (look at curves labeled by '—1.25' and '—1.5', respectively). This really surprising observation 
indicates that the Hall-MHD kink surface waves for that value of rj (= 4) are unstable only for negative relative 
Alfvenic Mach numbers A/a ^ — 1 - their normalized growth rates are plotted in Fig. [T3] This example shows us 
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that we must be very cautious with stating some very general criteria for the Kelvin-Hclmholtz instability onset with 
Hall-MHD waves propagating in flowing solar flux-tube plasmas. 

4 Conclusion and outlook 

Let us now summarize the basic results obtained in this study. In investigating the wave propagation along a jet 
moving with respect to the environment with a constant speed U we had to take into account the influence of two 
factors: (i) the Hall term in the generalized Ohm's law, and (ii) the flow itself. The combining effect of these two 
factors can be expressed as follows: 

— The Hall term generally limits the range of propagation of the wave modes not only for static tubes/layers but 
also for jets with positive Alfvcnic Mach numbers. The limiting normalized wave number, Kn^^it, is specifled by 
two plasma parameters: the densities ratio of the two plasma media (outside and inside the jet), 77, and the Hall 
parameter, e [look at Eq. (jlSp J. If this is the exact value for the waves propagating on a static flux tube, it is 
approximately the same for the waves travelling along a jet. One should emphasize that all dispersion curves, at 
relatively large Ks, lie on the left to the dispersion curve corresponding to Ma ~ (see Figs. [T] and [6]). However, 
when Ma becomes negative the real part of the phase velocity of the eigenmodes (kink and sausage waves) is forced 
to go beyond that limiting wave number in a region where the wave becomes unstable (or if you prefer, overstable). 
The instability which occurs should be of Kclvin-Helmholtz type. It is rather surprising that the instability onset 
starts with relatively large growth rates gradually decreasing with increasing the modulus of the Alfvenic Mach 
number (look at Figs. [51 and [TU)) . We note that such growth rates in the extended range of the waves' propagation 
was obtained for the same plasma-jet configuration for a different value of the parameter rj {— 0.64) [50| . It is 
worth mentioning that for negative Alfvenic Mach numbers we actually have (for each mode) two dispersion curves 
merging at A'umit, but unstable for negative Ma are the dispersion curves situated on the right to the cusp (see 
Fig. [3] and [5]). Our conclusion that the Kelvin-Helmholtz instability onset is only possible for negative relative 
Alfvenic Mach numbers is in agreement with a similar inference of Andries and Goossens . 

— It seems that the critical relative Alfvenic Mach number at which the Kelvin-Helmholtz instability starts, given 
by Eq. (|17p . is unfortunately not applicable, as a rough estimation, in the Hall magnetohydrodynamics. Notwith- 
standing, it is still useful because yields that value of the negative Alfvenic Mach number which is associated with 
a dramatic change in the shape of the waves' dispersion curves. 
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— A safely general conclusion is that the HaU current keeps stable the surface modes travelling in flowing solar plasmas 
within the dimensionless wave number range between and /iumit for each relative Alfvenic Mach number. An 
instability of the Kelvin-Helmholtz type is possible only at negative Alfvenic Mach numbers in a wave number 
range lying beyond the i^nmit and it (the instability) starts at some critical normalized wave number depending 
on the magnitude of A/a- The maximum instability growth rate is largest for small in magnitude Alfvenic Mach 
numbers gradually decreasing with the increase of |A/a|. The instability can turn off as |Ma| reaches some value 
depending on the magnitude of the parameter t] (look at Figs. [T2l and [T3ll . 

This study can be extended in two directions. The first one is to consider finite-valued sound speeds. In that case, 
however, the waves' dispersion equations become rather complicated and they can be solved (looking for complex 
roots) only numerically which is a highly difficult task. Any way, one can expect that the plasma compressibility will 
not substantially change the Kelvin-Helmholtz instability pictures. Moreover, one can state that the sausage and kink 
waves represented in Figs. 5 and 7 in Ref. [1^ are definitely stable because with 77 w 0.6 and e ~ 0.4 the value of 
^limit, beyond which one can expect the instability onset, is 3.16, while the waves' propagation in that paper was 
examined till kxQ = 2.5 only. The second direction is to study a more realistic geometry, for instance, cylindrical one, 
and conduct the investigations with appropriate observable plasma and magnetic field parameters. This is in progress 
and will be reported elsewhere. 

We do believe that the present results might be useful in studying wave turbulence in the solar wind as well as in 
solving other problems associated with wave propagation in structured/spatially bounded magnetized plasmas. 
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Fig. 1. (Online colour) Dispersion curves of kink Hall-MHD waves travelling along an incompressible flowing plasma layer for 
positive values of the relative Alfvenic Mach number Ma - all waves are stable. For details in curves' labeling see Fig. |2l 
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Fig. 3. (Online colour) Dispersion curves of kink Hall-MHD waves travelling along an incompressible flowing plasma layer for 
negative values of the relative Alfvenic Mach number Ma- All the curves lying on the left side of the vertical merging line/cusp 
represent stable waves' propagation. The waves become unstable only on the family of curves labeled by '—0.5' and ' — 1.45'. For 
details in curves' labeling see Fig. U 
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Fig. 4. (Online colour) A zoom of the bottom part of Fig. |3] The dispersion curves situated between the curves labeled by 
' — 1.25' and ' — 1.45' correspond to A'Ia equal to —1.3, —1.35, and —1.4, respectively. 
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Fig. 5. (Online colour) Growth rates of unstable kink Hall-MHD waves travelling along an incompressible flowing plasma layer 
for negative values of the relative Alfvenic Mach number Ma in the short- wavelength region (beyond the iCiimit). 
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Fig. 6. (Online colour) Dispersion curves of sausage Hall-MHD waves travelling along an incompressible flowing plasma layer 
for positive values of the relative Alfvenic Mach number Ma ~ all waves are stable. For details in curves' labeling see Fig. [T] 



22 



1. Zhelyazkov: Hall-MHD waves and instabilities in flowing solar plasmas 




Fig. 7. (Online colour) A zoom of the bottom part of Fig. |6l The dispersion curve corresponding to Ma ~ 1.315 divides the 
rest curves into two types, notably "simple" dispersion curves (for Ma < 1.315) and much complex families of dispersion curves 
(for Ma > 1.315). 

2.5 



^ 1.5 
DC 



0.5 



1 


1 


1 1 










|kl.45 


-■\ 




-1.3 






1 


1 1 


-0.5^'~^^ 



10 



kx 







Fig. 8. (Online colour) Dispersion curves of sausage Hall-MHD waves travelling along an incompressible flowing plasma layer 
for negative values of the relative Alfvenic Mach number Ma- All the curves lying on the left side of the vertical merging 
line/cusp represent stable waves' propagation. The waves become unstable only on the family of curves labeled by '—0.5' and 
' — 1.45'. For details in curves' labeling see Fig. [O] 
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Fig. 9. (Online colour) A zoom of the bottom part of Fig. [H The dispersion curves sandwiched between the curves labeled by 
' — 1.3' and ' — 1.45' correspond to Ma equal to —1.35 and —1.4, respectively. The curves situated at the left bottom corner of 
the dispersion diagram represent backward stable sausage surface waves. 
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Fig. 10. (Online colour) Growth rates of unstable sausage Hall-MHD waves travelling along an incompressible flowing plasma 
layer for negative values of the relative Alfvenic Mach number Ma in the short- wavelength region (beyond the -/^umit). 
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Fig. 11. (Online colour) Dependence of the normalized instability growth rate on the ratio po/po = at a fixed dimensionless 
wave number {kxo = 8.5) for various relative Alfvenic Mach numbers in the range between —1.45 and —0.5 The beginning 
of the horizontal axis starts at = 0.2. The family of curves sandwiched between the curves labeled by ' — 1.25' and ' — 1.45' 
correspond to relative Alfvenic Mach numbers equal to —1.3, —1.35, and —1.4, respectively. 
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Fig. 12. (Online colour) Dispersion curves of kink Hall-MHD waves travelling along an incompressible flowing plasma layer for 
negative values of the relative Alfvenic Mach number Ma and rj — 4. All the curves lying on the left side of the vertical merging 
line/cusp represent stable waves' propagation. The waves become unstable only on the family of curves labeled by '—0.25' and 
'— 1'. The curves with labels '—1.25' and ' — 1.5' represent stable waves' propagation. 



